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Abstract. The formulae for calculating jet fragmentation momentum, (jf'), and 
parton transverse momentum, (fc|,), and conditional yield are discussed in two particle 
correlation framework. Additional corrections are derived to account for the limited 
detector acceptance and inefficiency, for cases when the event mixing technique is used. 
The validity of our approach is confirmed with Monte-carlo simulation. 



1. Introduction 

In pp, pA or AA collisions at RHIC, the high pt part of the hadron spectra is dominated 
by the two-body hard-scattering process. In this process, the two scattered partons 
typically appear as a pair of almost back-to-back jets. The properties of the di-jet system 
can be characterized by the jet fragmentation momentum, jrp, the parton transverse 
momentum, kx, and the fragmentation function, D{z). Specifically, jy, kf and D{z) 
describe the spread of hadrons around the jet axis, the relative orentation of the back- 
to-back jets and the jet multiplicity, respectively. 

Traditionally, energetic jets were reconstructed directly using standard jet finding 
algorithms m 12]. In heavy-ion collisions, due to the large amount of soft background, 
direct jet reconstruction is difficult. Even in pA or pp collisions, the range of energy 
accessible to direct jet reconstruction is probably hmited to pr > 10 GeV jB], below 
which the 'underlying event' background contamination become important. The 
situation is even more complicated for finite acceptance detectors like PHENIX due 
to the leakage of the jet cone outside the acceptance. 

The two-particle-correlation technique provides an alternative way to access the 
properties of the jet. It is based on the fact that the fragments are tightly correlated 
in azimuth (j) and pseudo-rapidity t] and the jet signal manifests itself as a narrow peak 
in A0 and Arj space. The jet properties are extracted statistically by accumulating 
many events. This method was initially used in 70's in searching for jet signals in pp 
collisions at CERN ISR [H El Ej- It overcomes problems with the underlying event 
background, probes the jet signal at lower pxi and has recently excited renewed interest 
at RHIC 0012]. 
§ jjia@shang.nevis.columbia.edu 



Probing jet properties via two particle correlation method 



2 



Experimentally, jet measurement is challenging due to detector inefficiency and 
limited experimental acceptance. A certain fraction of jet pairs is lost either because 
the track is not found or because part of the jet cone falls outside the acceptance. The 
average jet pair detection efficiency can be estimated statistically using event mixing 
technique. The main goal of the paper is to establish the procedures for measuring the 
jet shape and jet multiplicity independent of the detector efficiency and acceptance. 

The discussion is split into three sections. In Section|21 we lay out the formulae for 
extracting Jt, kx and conditional yield. In the interest of space, the reader is encouraged 
to see Ref. CHI for more details on the definition of the variables. In SectionI2l we 
discuss the event mixing technique necessary in correcting for limited acceptance and 
inefficiency, where we use the PHENIX detector as an example. We discuss separately 
the two dimensional [Ar], A0) and one dimensional (A0) correlation and derive the 
normalization factors for the conditional yield in each case. In SectionlUwe verify the 
procedure of correcting for finite acceptance with Pythia simulation. 



2. Some formulae 



2.1. Formulae for Jt , kr 

FigHillustrates the single jet fragmentation (left panel) and away side jet fragmentation 
(right panel). The jet fragmentation momentum, Jt, and parton initial momentum, kr, 
determine the relative orientation of the fragmented hadrons. Jt and kx are vectors, 
their projection in the azimuthal plane perpendicular to the jet direction are denoted 
as and kxy- For single jet fragmentation, if we denote Acp, (ptj, and (paj as the 
angles between trigger-associated, trigger-jet and associated-jet, respectively, then the 
following relations are true: 

A0 = (ptj + (paj , sin(A0) = , 

PT,asso 

. / , N jTy ,trig _ ■ r i \ 3Ty ,asso \ 

Sm{(ptj) = = Xj^trig , Sm{(f)aj) = — = Xj^asso ■ (1) 

PT,trig PT,asso 

Here Pout,N is the component of associated particle px {pT,asso) perpendicular to the 
trigger particle pr {pT,trig)- Assuming (ptj and (paj are statistically independent, we have 
(cross terms average to 0), 

^sin^ A0^ = ^sin^ (ptj cos^ 0aj) + (sin^ (paj cos^ (ptj'^ (2) 

Similarly for the far-side correlation, we have 

A0 = (j)tj + (f)aj + 4>jj , sin(A0) = , (3) 

PT,asso 

where (pjj represents the angle between the two jets. Expanding sin^ A0 and dropping 
all cross terms (which average to 0), we get 

2 \ i\ 1 1 ■ X X X \^\ X 1 1 ■ X X X \^ 



<^sin A0y = (^{sin (j)tj cos (j)ajCos(f)jj) J + (^{sm(f)aj cos (f)tjCos(j)jj) J + 

(|(sin cos (j)ajCOS(f)tjf^ + {^(sin (ptj sin (pajSirKpjjf^ (4) 
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Let's define ZtHg = Pt, trig Ipt, jet and Xk = V2kTy/pT,jet = V^kTyZtrig/pT,trig- If 
is small, or prjet > ^t^, then sm{(f)jj) ^ {kTy,trig + ^T^^asso) /PTjet, ^-e. /ct only affects 
the jet direction. Thus statistically we obtain 

(.mU.)>-2(-^\ = (4 . (5) 

\ PT,jet I 

Substituting the sin and cos terms from Eq. ^ and El into Eq. |21and Eq. EJ we obtain the 



equations for the RMS value of jTy and kTy (for a given variable x, {x)rms = y i^^))^ 



TyZtrig 



RMS 



Plut^p) - {plut,N) (1 - (4)) 



\ 2 ((^^i) i^^'^Xtrig) {^^'j,trig} + ^ (x^^f-^.^ 



(7) 

"],trig/ 

where Xh = PT,asso/pT,trig- Assumiug the A0 follows gaussian statistics, a simple Taylor 
expansion connects pout with the jet width, a: 

{pint) = {PT,assoSin^^^) 

Pt,„.„) l^i""' - -^1 (8) 
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By replacing the Pout terms in EqO and Eq01 we can derive the relation between jet 
width and RMS value of and . 

Since Eq. IHland Eq. [3 contain variables Xj^trig and Xk that depend on and kxy, 
we have to calculate {jTy)RMS and {kTy)RMS iteratively. However when trigger and 
associated particle pr are much larger than typical jt value, the near side jet width 
is small and Xj^trig ~ 0. Hence Eq. El can be simplified to, 



\^^y)RMs-\ l + {xl) ^' 

For the far-side correlation, if the trigger pt is much larger than kj^y and jxy, then 
{xD ^ 0, (x'j) ^ and Eq. [7| reduces to 



Pout,F) \Pout,N 



\ 2 {xD 



(10) 



2(4) 



\l {PT,assocY Sin^ aF-{l+ {xD) (jTy) 



2 

RMS 



Previously, Ref jU] has derived the approximate relations between the j^, kf and 
the measured jet width in the near side (div) and the away side (cr^) ((|a|) is the mean 
of the |a|, {\a\f = 2/Tx{a^)): 

(lirj) = ^{pT,asso) <TN j ^ I + {xhf (11) 

{\kTyZtrig\) = ^1^^^ ^ lj)T,asso)'^ ]J^CrA - (1 + (Xhf) (brj)^ (12) 

One can see that Eq^^and EqJ12lare very close to our approximation EqEland EqllOl 
To quantify the difference between EqlHlIZIand Eq JTTlIT^ we performed the correlation 
analysis with vr^ — correlation from the Pythia event generator ^T] and compared 
the Jt kx value calcualted from the two sets of formulae. The results are shown in 
Fig|21 there are good agreements for jxy, except at low pT,asso where the become big 
and EqlHlhas to be used. On the other hand, the (ztrigkr^ value from this analysis is 
10% lower than the one calculated by Eq|7| The difference could be caused by small 
angle approximation used in EqEl but may also due to the difference between RMS and 



mean: for given variable f in a finite range, y (v"^) > (v). There is a dropping trend in 

the (ztrigkT^ value as function oipx^asso in both cases. This trend can be explained by 
the trigger bias effects ^2j, which we briefiy touch on in the next section. 



2.2. Formula for yield 

The jet fragmentation function, D{z) or D{pt), represents the associated yield per jet, 
D{z) = l/Nj,tdNasso/dz, or D{j>t) = l/NjetdNasso/dpr (13) 
where D{pt) = D{z)/pT,jet, when the jet energy prjet is fixed. 
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Figure 2. The comparison of the two sets of formulae for Jt^ RMS value (left panel) and Ztrigkr^ 
RMS value (right panel). Values are calculated from jet width given by Pythia simulation. 



If the two particles come from the same jet, the two particle multiplicity distribution 
can be described by the di-hadron fragmentation function, 

D{zr,Z2) = l/N^,td^Nh/dz^dz2 . (14) 

The conditional fragmentation function, where the pt of one particle is fixed, can be 
expressed by, 

D{zi,Z2) ^ 1 d-'Nh I 1 dN^ 

D{zi) NjetdZidZ2/ Njet dzi 

Experimentally, we measure jet properties from two particle azimuthal correlation. 
The jet signal typically appears as two distinct peaks in A0, i.e. at A0 = for same 
jet and at A0 = vr for away side jet. By fitting the peaks with gaussian function, we 
extract the associated hadron yield per trigger or the conditional yield {CY), 

1 d'^N 

Cr(A0, Ar^) = ^^^^^ (16) 
^ " Ntrig dA(t)dAr] ^ ' 

Where 'Ntrig is the number of trigger particles, and ^/^^^ is the number of jet pairs 
in a fixed A0 and Ar] bin. Since Ntrig oc D[ztrig) and ^^^^^ oc D^ztrig, Zasso ), the CY 
integrated over the same jet (around A0 = 0) directly corresponds to the conditional 
fragmentation function. Unfortunately, the fragmentation variable z is not directly 
accessible in the correlation method. Instead, CY is often expressed as function of a 
different variable xe defined as xe =Pt ■ PT,trig/\PT,trig\'^ = Zasso/ Ztrig- 

The di-hadron fragmentation function (Eq lT^ . conditional fragmentation function 
(EqlT^ and conditional yield (Eqll6p are defined for the near side correlation for 
which the two particles belong to the same jet. They can be extended to describe the 
correlation of the two particles from the back-to-back jets. In this case, Eq^J represents 
the hadron-triggered fragmentation function similar to the one used in ^31- In a naive 
parton-parton scattering picture and for fixed jet energy, the fragmentation of the away 
side jet and triggering jet can be factorized, i.e. D{zi,Z2) = '^^l"^ ivH '^^dll"" ■ Thus 
the CY is related to the unbiased fragmentation function by a scale factor, (ztrig), 



1 dN, 



CY{xE) = ^^ = Zt„gDiz) . (17) 

Ntrig dXE 



Probing jet properties via two particle correlation method 



6 



(ztrig) was found to be around 0.75-0.95 at ISR |Sj for high pt leading pions (4-12 GeV/c) 



In reahty, due to the intrinsic kr or radiative corrections, the independent 
fragmentation assumption for parton-parton scattering is not strictly valid. In addition, 
for given given pT,trig, the original jet energy is not fixed, but depends on the PT,asso, 
which implies that Ztrig also depends on pT,asso- This can be easily understood from the 
fact that the fractional momentum can't exceed 1: < 2; < 1, while xe is not bounded 
from above: < x^; < 00. Simple simulation indicates that Zt^ig is relatively stable when 
PT,asso < PT,trig, but quickly dccrcasc when pT,asso is larger than pT,trig [12] • In summary, 
these bias effects, caused by the requirement of a trigger hadron on the opposite side, 
makes EqEI at best an approximation. The trigger bias effect is also responsible for 
the decreasing of (ztrigkTy) seen in Fig|2l 



3. Corrections for limited acceptance and efficiency 
3.1. CY in A(j) and Ar] 

For a detector with limited aperture, a fraction of the jet cone falls outside the 
acceptance. The fractional loss depends on the jet direction: the loss for jet pointing 
to the corner of the detector is larger than jet pointing to the center of the detector. 
Assuming the jet production rate is uniform in azimuth direction (p and pseudo-rapidity 
rj, we can construct an average pair acceptance function using event mixing technique 
(One trigger particle is randomly combined with an associated particle from a different 
event.) 



mo represents the background level when the associated particle is not constrained. 
Given the limited PHENIX rj coverage, it is safe to assume that mo is constant |14j . 
Acc{A(j), Arj) is the pair acceptance function, which represents the probability of 
detecting the associated particle when the trigger particle is detected. 

The PHENIX ideal single particle acceptance is shown in FigO^ : G 
[-33.75°, 56.25°], [123.75°, 213.75°] and \ri\ < 0.35. The tracking efficiency has been 
assumed to be 100%. The corresponding pair acceptance can be constructed by 
convoluting the single acceptance for trigger and associated particle and is shown in 
FigOb for Arj and FiglHt for A0. Interestingly, although the PHENIX detector has only 
TT coverage in azimuth, the pair acceptance actually is sensitive to the full range in A0, 
^A(/. = 277; similarly the pair coverage in Arj is also doubled relative to single particle 
rj acceptance, Q^ri = 1-4. The pair phase space is 27r x 1.4 = 2.87r, four times that 
for single particle acceptance, vr x 0.7 = 0.77r. However, this increase is compensated 
by the decrease in overall pair efficiency, which is only 25%. There is only one point 
(A0 = 0, Arj = 0) where Acc = 1 (the probability of detecting the associated particle 
when the trigger particle is detected is 1). 



and scales with xt in y/s [S]. 




(18) 
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Figure 3. The PHENIX ideal 
single particle acceptance (top left) 
and corresponding pair acceptance in 
A77 (top right) and A(p (bottom). 



The foreground distribution is modulated by the same pair acceptance function 

Acc{A(j), Ar]){Xmo + jet{A(f), Ar])) 



dA(j)dAr] 



X- 



+ 



(19) 

'dA(j)dAT] ' dA(j)dAr]' ^'^^^ 
By dividing the foreground by the mix distribution, the Acc cancels out, and we 
are left with a constant background plus the jet signal. This ratio has correct jet shape, 
but the magnitude is off by factor of l/mo- Let's assume the original number of trigger 
{\r]\ < 0.35) and associated particles (within |A?7| < 0.7) per event are Utrig and Uasso, 
and those within PHENIX acceptance per event are nftrig and nfasso, respectively. Then 
the sum rule for original and acceptance filtered mix event pair distribution are. 



j2 ]\Tmix 

dA4>dA^^£^ 
nA4>x^Ar, dAcpdAr] 

dA<f)dA7] 
^A^x^A,, dAcpdAr] 



^events'^trig'^asso 



events^^^ trigf^^ asso 



(21) 



(22) 



where Ntrig 



Neventsntrig is the total number of triggers. From EqUHJES we obtain the 
following relation between the original (true) jet CY and the measured CY: 

1 d^N^o'' 1 1 d'^N^''/dA(j)dAr] 



Na dA4)dAr] Nu e« 

so/lT'ass 



J JsiA^xnAr, 



" /dA(f>dAri 



IT'^asso/ IT-asso ^ single / ^ Arf 

^single rcprescuts the single particle correction to 27r in and 1 unit in 77. 



(23) 
(24) 



3.2. CY in A(j) 



Building two dimensional correlation requires high statistics for event mixing. Instead, 
correlation function is often built as function of A(f) only, by integrating the 2D 
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correlation function over Arj. The following relations are true, 



1 dN^'' 1 f d^Nt 



Ntrig dAcj) Ntrig JnAr, dA(j)dAri 
1 dN^^^ 1 r d^iV^'^* 



dAT]——^ , (25) 



dAr]——— , (26) 



Nftrig dAcj) Nitrig J^Ar, dA(j)dAT] 

dA(j) Jqav ^ dAcpdArj ^ ^ 

To relate the measured A0 distribution fEgOUj) with the true A0 distribution 
(EqEni), we require the following two assumptions: 

• The jet signal can be factorized in Ar] and A0, i.e. jet{Ari, A0) = gi{A(j)) x (72 (Ar^) 
Typically, the jet signal can be approximated by, 

jet{A(P,Ari) = {C^e ^ + C^e ''"I ) x (Cae ^ + J{Ai])) .(28) 

In pp or pA collisions, near side jet widths in A0 and Ar] are the same, a\ = a^- 
On the away side the jet distribution in Ar] {J{Ari)) is typically very wide due to 
the difference in momentum fraction, x, of the two initial partons. 

• The pair acceptance function can be factorized in A0 and Ar], i.e Acc = Acci{A(f)) x 
Acc2{At]). This condition is satisfied if the single particle efficiency factorize in (p 
and 77. 

with these two assumptions, using Eq|2Sl we can derive following relation that connects 
EqESland EqEEl 

1 dNt^ 1 dN^^'/dA(j) 



Ntrig dA<P N/trigeasso ^A^dN'^^/dA^ 

L dA(j>dN-^^^ /dA<j) 



(29) 



_ InA,d^vAcc2{A7])g2{Ar]) 

^single /nri\ 

^ dAr]Acc,{Ar]) J^^^ dA^g^iAr]) " ^ > 

Comparing with Eql24| we sees that EqlHIl have two more terms that are related to the 
Ar] correction. The first term is the pair acceptance weighted average of jet signal over 
Hat]- The second term (/q^^ dAr]g2{Ar])) corresponds to the fact that only a certain 
fraction of jet signal falls in the l^Ary- 

EqI2ni and EqI2ni are equivalent: both correct the jet yield to the true yield in 
^A??, i-e \Ar]\ < 0.7. To account for the loss of jet yield outside ^Ar?, an additional 
extrapolation factor is needed. Because the jet shape in Ar] are quite different for the 
same side and away side, we use two different approaches. For the same side we assume 
the jet width are equal for Ar] and A0, i.e. o"i = (72 [HI; and simply extrapolate g2 
(guassian) to the full jet yield. On the away side, since the jet width in Ar] is much 
broader than the typical PHENIX pair acceptance, we assume g2 is constant in f^Ar; and 
no additional extrapolation is required. Under such approach. Eg ]30l becomes. 
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dAr}Acc2{Ar]) 



full jet yield at the same side 



(31) 



asso 



yield in fi^r^ at the away side 



Since Acc2 typically has a triangular shape, this correction is easy to evaluate. 
4. Pythia Simulation 

The extraction procedures for jet width and CY are verified with Pythia event 
generator ^T], using ir^ — h^ correlation. A single particle acceptance filter is imposed to 
randomly accept charged particles according to the detector efficiency. FiglH shows the 
PHENIX style two dimensional single particle acceptance filter used in the simulation. 
The average efficiency in 27i in azimuth and 1 unit of pseudo-rapidity is egingie = 0.256. 



We generated 1 million Pythia events, each required to have at least one > 6 
GeV/c charged pions. To speed up the event generation, a cut of > lOOGeV"^ on 
the underlying parton-parton scattering is required. These events were filtered through 
the single acceptance filter. As an approximation, we ignore the pt dependence of 
acceptance. The same event and mixed pair A0 distributions were then built by 
combining the accepted vr^ and charged hadrons, where the trigger vr^ is selected to 
be 6 < PT,trig < 10 GcV/c. The jet width and raw yield were extracted by fitting the 
with a constant plus double gaussian function. The raw yields were then 
corrected via Eql^to full jet yield for same side and to true yield in | At^I < 0.7 for away 
side. Meanwhile, we also extract the true CY and jet width without the acceptance 
requirement. The comparison of the CY and jet width with and without the acceptance 
requirement are shown in Fig|S| In the near side, the corrected yield (top left panel) 
and width (bottom left panel) are compared with those extracted without acceptance 
ffiter. In the away side, the yield corrected back to lAr^l < 0.7 (top right panel) and 




-1 



^ (rad) 



Figure 4. A typical 
PHENIX single parti- 
cle acceptance/efficiency 
map for charged hadrons 
that was used in the sim- 
ulation and plotted as 
function of and Arj. 
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width (bottom right panel) are compared with those extracted without acceptance filter. 
The data requiring acceptance filter are always indicated by the filled circles while the 
expected yield or width are indicated with open circles. 
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- ture yield 
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- ture yield in |Ar||<0.7 
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Figure 5. The comparison of near 
side yield (top left panel), near 
side width (bottom left panel), 
away side yield(top right panel) 
and away side width (bottom 
right panel) as function of pT 
of charged hadrons. These are 
obtained for tt^ — correlation 
from Pythia, with trigger pion 
from 6-10 GeV/c. The open circles 
represent the quantities calculated 
with the acceptance filter show in 
FigH 



The agreement between the two data sets can be better seen by plotting the ratios, 
which are shown in FigEl The yields agree within 10% and the widths agree within 
5%. Since (j^ ), (k^ ) are derived from the jet width(Eq|21E)), the agreement in width 



naturally leads to the agreement in the \3Ty)i \^Ty)- notice that there are some 

systematic difference in the comparison of the yield at low PT,asso- This might indicate 
that the gaussian assumption is not good enough when the jet width is wide and the 
extrapolation for IAt^I > 0.7 become sizeable (At PT,asso = 0.5 GeV/c, the jet width 
(Tat = 0.5(rad), and the extrapolation is about 20%.). 
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Figure 6. The ratio of the 
jet width or corrected yield 
obtained using event mixing 
method to those obtained 
without acceptance filter. 



5. Conclusion 

The formulae on Jt, kx and CY are discussed in two particle correlation framework. 
A more general definition of {jTy)RMS is found to be consistent with previous 
approximation, but the {kTy)RMS is lower by up to 10%. We have also demonstrated 
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that the event mixing technique can reproduce the jet width (thus the Jt and k^) for 
a hmited acceptance detector. With a correction factor that takes into account the 
hmited Arj acceptance, we can also reproduce the CY. Based on a Pythia simulation 
in PHENIX acceptance, our procedure can reproduce the CY within 10% and the jet 
width within 5% at both the same side and away side. 
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